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Moment densities of super-Brownian motion, and a 
Harnack estimate for a class of X-harmonic functions 



Thomas S. Salisbury and A. Deniz Sezer ^ 

York University and University of Calgary 



Abstract: This paper features a comparison inequality for the densities of the moment measures of 
super-Brownian motion. These densities are defined recursively for each n > 1 in terms of the Poisson 
and Green's kernels, hence can be analyzed using the techniques of classical potential theory. When 
n = 1, the moment density is equal to the Poisson kernel, and the comparison is simply the classical 
(-H ' inequality of Harnack. For n > 1 we find that the constant in the comparison inequality grows at most 

_ exponentially with n. We apply this to a class of X-harmonic functions H'^ of super-Brownian motion, 

introduced by Dynkin. Wc show that for a.e. in this class, H" {^) < oo for every ^. 



1. A moment density inequality 



. Let D he a bounded and smooth domain in W^^ d> 2 and let A = Yli=i '§^- Let goix, y) and ko{x, y) be 

respectively the Green's and Poisson kernels of D for the operator A. That is, for x E D, goix, ■) : D — {x} — > 
^ ■ [0, oo) is the unique C^{D - {x}) nC^{D- {x}) solution of 



Aw — 26^ on D 
u — on dD. 

, Then for z £ dD, ko{x, y) € [0, oo) is the derivative of goix, z), considered as a function of z, in the direction 

■ of the inward conormal Uy to dD at the point y € dD. For fixed y € dD, the function u — k£i{-,y) satisfies 

. Au = VI on D 

^ ; M on - {y}. 

■ For y £ D, linix^y gD{x,y) = oo. Similarly, for fixed y g dD, linix^y kD{x,y) = oo, provided the limit is 

■ taken non-tangentially. (See Doob (1984) for these and other basic results in potential theory.) The notation 
D' D will mean that D' has compact closure which is contained in D. 

Let ri > 1, and let zi,...,Zn be distinct points on dD. We will define a function pd{x, zi, . . . , z„), 
X <E D recursively, by defining a family of functions pz)(x,za), where A — {ii, . . . ,ik} C {l,...,n} and 
ZA = (zji, ■ • ■ , ZjJ- Set 



and for \A\ > 1 set 



PDix,ZA) = kD{x,Zi) if A = {i}, (1) 



PDix,ZA) = y] / 9D{x,y)pD{y,ZB)pD{y,ZA-B)dy. (2) 



BcA.B^Hl.A 



We let pd(x, zi, . . . , Zn) = p{x, za) for A = {1, . . . , n\. It is known that pd{x, zi, . . . , z„) is a finite valued 
function oi x G D. See, for example Dynkin (2004a) or Salisbury and Verzani (1999). In section 3 we will 
interpret p as a moment density. Unless needed for clarity we will drop the subscript D in pu, k]j and go. 
The main estimate of this section is the following theorem: 



^Both authors are supported in part by NSERC. This version June 29, 2012. MSG 60J68 (primary); 60J45, 31B05 (secondary) 
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Theorem 1. Let D be a smooth and bounded domain in W^, d> 2. For any compact C C D, there exists a 
A > depending only on C and D such that 

p{x, Zi,...,Zn) < A"'p(xo, Zi,..., Zn) 

for all X and xq Cz C and for all distinct zi, . . . , z„ G (9Z3". 
Wc will describe the implications of this result in section 2. 

One should compare the above estimate to Harnack's inequality: Indeed, when n = 1, p{x,z) = k(x,z). 
The theorem in this case therefore follows from Harnack's inequality, which asserts that there exists 9c 
depending only on C and D s.t k{x, z) < 9ck{xo, z) for all x, xq G C and z G dD. 

The proof for general n will be given in Section 4. To see where the main challenge in the proof arises, 
consider the case n = 2. Then, p(x, zi, Z2) = 2 g{x, y)k(y, zi)k{y, z^) dy. Due to Harnack's inequality, there 
exists a 6d' such that g{x, y) < 9uig{xo, y) for any y G D — D', where C C D' , and D' d D. If this inequality 
were to hold on the entire set D it would be easy to compare p{x,zi,Z2) with p(a;o, zi, Z2). However this 
is not the case, as 17(3;, y) blows up when y is near x. Since the integral defining p is taken over the entire 
domain D, in order to compare p{x, •) with p{xo, ■) we need finer estimates for g{x, y) when y is near x, and 
for fc(y, zi)k{y, Z2) when y is near zi or Z2. 

The proof follows an inductive argument. We will obtain compact subsets Co D Ci D . . . D C of D. Then 
for each n, we will prove that 

p(a;, zi,...,z„) < A"p(a;o,zi, . . . ,z„) 

for all x^xq G Cm assuming 

p{x, zi, . . . , z„_i) < A""V(2;o, ^1, • ■ • , Zn-l) 
for all x,xo G Cn~i- We do this by decomposing Jj^ g{x,y)p{y, ZB)p{y, ZB<')dy as 



9{x,y)p{y,ZB)p{y,ZB'-)dy + I g{x,y)p{y,ZB)p{y,ZB'-)dy 

B{x,S„) Jd-B{x,S„) 

where (5„ = jf^- Here S is the diameter of the domain and iV is a certain constant strictly greater than 1. 
For the first term, we use the induction hypothesis and the fact that g{x,y) behaves like \x — y\'^~'^ when y 
is near x to show that 



/ 9{x,y)p{y,ZB)p{y,ZB'=)dy < ^ [ 

Jb{x,S„) ^ Jb{xo,5'„[ 



ai^o, y)p{y, ZB)p{y, ZB'')dy 



for some A > 0. For the second term we use the 3-G inequality. In the case d > 3 this gives a 9 depending 
only on D such that 

9{x,y) ^ - yl^-'^ + \x - xo\^-'^ 



yi^o^y) gix,XQ) 

The critical factor on the right side is |a; — yp^'^, which results in a bound for gl^^^^y^ on D — B{x, Sn) of the 
order of n'^'^ which can in turn be bounded by A" if A is large enough. This implies that 

f A" f 

/ 9{x,y)p{y,ZB)p{y,ZB-^)dy < — g{x,y)p{y,ZB)p{y,ZB'^)dy. 

JD-B(x,5„) ^ Jd-B{x,5„) 

Putting these two pieces together will prove the theorem. 

The paper is organized as follows. In Section 2 we will describe an application of Theorem 1 to a finite- 
ness question for an interesting class of extended X-harmonic functions related to super-Brownian motion 
(originally introduced by Dynkin). Section 3 contains the proof of this result, and the proof of Theorem 1 is 
given in Section 4. 
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2. X-harmonic functions 



The functions pd{x, zi, . . . z„) will arise as moment densities of super-Brownian motion. More precisely, we 
consider a super-Brownian motion on some domain E (larger than all domains D of interest), as represented 
by Dynkin in terms of exit measures. See e.g. Dynkin (2002) for the definition and construction. This is a 
branching-exit-Markov system {Xo^ P^j) indexed by sub-domains D of E, and by finite measures ^ supported 
in E. Under the probability law P^, each Xd is a random finite measure, referred to as the exit measure from 
D of super-Brownian motion started with measure ^. Heuristically, in a branching particle approximation 
to super-Brownian motion, it would be approximated by a multiple of the empirical measure of the particles 
frozen when they exit D. More formally, it has the following basic properties: 

(a) Exit property: P^{Xo{E>) = 0) = 1 for every /i, and if n{D) = then P^{X]j — fi) ~ 1. 

(b) Markov property: If F > is measurable with respect to the cr-algebra J-cd generated by Xw, D' C D 
and Z > is measurable with respect to the cr-algebra J^z)D generated by Xd", D" D D then 

P^{YZ) = P^{YPx^Z). 

(c) Branching property: For any non-negative Borel /, Pfj,{e^^'^"'^'^) = e^^^^''^"^'^ where 

VDf{y) = -logPy{e-^''"^f^) 

and Py = Ps^ . 

(d) Integral equation for the log-Laplace functional: Vnf solves the integral equation 

where Go and Ko are respectively Green and Poisson operators for Brownian motion in D. 

The operators referred to above are defined as follows: If is a Brownian motion starting from x, under a 
probability measure n^;, then KDf{x) — ^xf{£,TD)^ where tu is the exit time from D. Likewise, Gnf{x) = 
nx(/g " fiS,t)dt). In terms of the kernels go and ko defined earlier, Gofix) = goix, y) dy and Kof{x) = 
J f{y)kD{x,y)a{dy), where cr is the surface measure on dD. 

Under certain regularity conditions on D and / (see e.g. Dynkin (2002)), the integral equation in (d) is 
equivalent to the boundary value problem 

^Am = 2u^ on D 
2 

u = f on dD. 

A recent research program, initiated by Dynkin, aims to build a Martin boundary theory for the above 
semi-linear pde employing ideas from the probabilistic construction of Martin boundary for the Laplace 
equation jAu = 0. In this new context, the analogue of harmonic functions are X-harmonic functions, 
defined as follows. Let Ai'jj denote the space of finite measures supported on a compact subset of D. Consider 
a function H : Ai'jj — > [0,oo]. H is called extended X -harmonic in D if 

P^{H{Xn')) = Hifi) 

for all D' D and /i s.t. /i is supported in D'. If, in addition H{pi) is finite for all jj, £ then H is called 

X -harmonic (this definition is due to Dynkin - see Dynkin (2004a)). 
Define 

P^,x„(^) := P^.{Xd eA.Xo^Q). 

Theorem 5.3.2 of Dynkin (2004b) shows that P^^.Xd ^^^d Pfj,2,XD ^-re mutually absolutely continuous. We fix 
xq G D, and take R = Ps^g,Xo a reference measure. Following Dynkin (2006b), we define 

ftp 

H^{ix) = HUt^) = ^^{v). (3) 
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Unless there is ambiguity about the domain D in question, we wih suppress the subscript D in H'^. 

Of course, for each fi this is only uniquely defined for R-a.e. v, so to proceed further we will need to 
specify more regular versions. The following is contained in Theorem 2 of Salisbury and Sezer (2012) (whose 
proof is a modification of that of Theorem 1.1 of Dynkin (2006b)), and establishes the existence of a jointly 
measurable version of H'^{n) that is extended X-harmonic for _R-almost all i'. We assume that D is a. bounded 
domain all of whose boundary points are regular. The smoothness assumption of Theorem 1 is not required 
here. 

Theorem 2. (Dynkin (2006b), Salisbury and Sezer (2012)) Let D be a bounded and regular domain in , 
d > 2. There exists a measurable function (fj,, v) H'^(fi) G [0, oo], defined for finite /i € M'fj and finite v 
supported on dD, and an R-null set Vg such that: 

(a) For each v ^ Yq, [i i— > H'^(fi) is extended X-harmonic. 

(b) For each /i, i— > H'^{ii) is a version of the density in (3). 

For each fixed /i, (b) implies that H'^{^) < oo for R-a.e. v. But the i?-null set where H'^ifJ-) = oo could in 
principle depend on yit, so one does not automatically obtain any v for which H^{^) < oo\/fi. In other words, 
it is not clear that any of the H'^ is X-harmonic. However, we can establish this using Theorem 1: 

Theorem 3. Let D be a bounded and regular domain in M^, d> 2. There exists an H'^{fi) as in Theorem 
2, and an R-null set Vi such that for v <^Vi, H'^ is X-harmonic, and < H'^{p) < oo for every fi. 

The above is the principle result of this section, and the main application of Theorem 1. We give the 
proof in section 3, and will spend the remainder of the current section setting the context for this result, and 
giving relations to other work. 

Note first that the statement of Theorem 1.1 of Dynkin (2006b) is stronger (ie. X-harmonic) than than 
what is verified in the proof there (ie. what we call extended X-harmonic). Obtaining the required finiteness 
estimate was the motivation for the current paper. 

If is a non-zero X-harmonic function one can proceed to define an H -transform of the law P^. Set 
Pf^ 0^) — H{fi) -Pfj-O^ H {-^ D' j) for any non-negative Y measurable with respect to J^u' = a{Xj=),D C D') 
where D' ^ D contains the support of fi. Because H is X-harmonic, the Markov property implies that these 
probability laws will be consistent as we vary D' , so can be uniquely extended to J-d- ~ <^{Xd' , L)' <s D), 
as in Section 1.3 of Dynkin (2006b). See also Theorem 2(d) of Salisbury and Sezer (2012). Note that the 
mutual absolute continuity of P^i,Xd and Pii2,Xo: for fii and /i2 compactly supported in D, implies that 
if H{iJ,) is non-zero from one ^ then it is non-zero for all fi. Thus the only obstacle to defining for H 
extended X-harmonic is the possibility that H^fi) = oo. 

Fix a E D, a ^ xq and let T-L"" be the convex set of all X-harmonic functions s.t. Ps^{H{X£)i)) = 1 for 
some D' with a £ D' D. The (minimal) Martin boundary is then defined as the set of extreme elements 
of "H". In Dynkin (2004a) it is shown that this set is independent of the reference point a. One motivation 
for introducing the Radon- Nikodym densities H^{ti) of Dynkin (2006b) was to try to recover the extreme 
X-harmonic functions by a deterministic limiting procedure from these densities, using the classical exit 
theory for Markov chains. What is actually known is a weak version of such a result. Specifically, Dynkin 
succeeded in showing that if H is an extreme X-harmonic function in D and Dn (s is a sequence of domains 
exhausting D, then iJ(/x) = lim„_j,oo Hd^" P^ almost surely. A stronger version is conjectured, namely 
that there exists a deterministic sequence (not depending on /i) and canonical versions of the ^^j^" (eg. 
having some type of regularity in Vn), such that H{^) = lim„_>oo H'^^{^) for every ^. Dynkin has pointed 

out that significant progress could be made if one knew that the densities -ff^"" were uniformly integrable. 
A subsequent paper Dynkin (2006a) focused on the densities if]^(/i) and the search for more analytically 
tractable representations. 

The densities H'j^ played a key role in Salisbury and Sezer (2012) as well, in particular, in the analysis of 
super-Brownian motion {Xd')d'ieD conditioned on its exit measure Xd- In addition to Theorem 2, that paper 
showed that P^f" is the conditional law of SBM given Xd = v, for R-a.e. v satisfying Hfy{fj,) < oo. Assuming 
that D is a regular domain, they express as an infinite sum of integrated polynomial-like functions of /x 
involving a family of fragmentation kernels if„(^, dvi, . . . ,Vn) and a family of potentials 7iy(-)i (see formula 
(6) below). This representation is then used to derive an infinite fragmentation system description for P^'^ . 
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In the current paper we will use the representation of Salisbury and Sezer (2012) to tie the densities 
H'^ to the functions p{x^ -z, ■ • ■ , Zn) and prove Theorem 3. This suggests that the comparison inequality of 
Theorem 1 may prove useful in investigating other regularity properties of H'^. In a subsequent paper, one 
of us (Sezer) will study the extremality properties of this family for smooth domains. As a future direction 
of research, we hope to use the comparison inequality to investigate uniform integrability properties of the 
sequences iJ]^" of Dynkin (2006b). 

We conjecture that a much stronger version of Theorem 3 is true. Namely that if H is any extended 
X-harmonic function, and H{^) < oo for some fi, then H{p) < oo for all fi. 



3. Proof of Theorem 3 



It is well known that Xo has an infinitely divisible distribution for each D. This property is tied to a 
construction of the probabilities in terms of another (now tr- finite) measure N^;, called the super- Brownian 
excursion law. Once again, there are finite exit measures Xd under Nx, for every D. Form a Poisson random 
measure 

i 

with intensity 7^p(A) = / Nx{Xd G A,Xd^ 0)^j.{dx). Then under P^, Xd has the law oij^x' = J X^{dx)- 
The n-dimcnsional moment measures of SBM are the following measures defined on ((?£>)": Let f{zi, . . . , z„) 
/i(zi) . . . fn{zn) where fi are positive Borel functions on dD, and define 

Nx,nMf) ^^x{{XdJi) ■ ■ ■ {XdJu)) (4) 
Pt^.nAf) = PMXdJi) ■ ■ ■ {Xd. /„)). (5) 

Theorem 4 (Theorem 5.3.1 of Dynkin (2004b)). Let x G D, and let /i be compactly supported in D. If D 
is smooth then both N^^n.D md p^^n.D have densities with respect to a{dzi) x a{dz2) x ••• x a[dzn) where 
<j{dz) denotes the surface measure on dD. Moreover, 

Nx,„,D{dzi, dzn) = p{x, zi, . . . , Zn)a{dzi)a{dz2) ■ ■ ■ cr(dz„) 

Pfj^^n.oidzi, . . .,dzn) = ^{n,p{-, zcj) ■ ■ ■ {p, p{-, zc^))a{dzi)a{dz2) ■ ■ ■ <y{dzn) 

7r(n) 

where irln) denotes the set of partitions {Ci, . . . , C^} of {1, ... ,n}. 

The proof of Theorem 3 involves only p^^n.D, not Nx,n,D- But the latter gives the natural context for 
directly interpreting p. 

Let V* be the space of nonzero and finite measures on dD. Let u{x) = N2;(X£) 7^ 0). For a regular domain, 
Salisbury and Sezer (2012) construct a family of functions {-ji, : D ^ (0, 00), v e y*} such that the mapping 
{y^ y) I— 7i/(i/) is measurable, each 7^^ is super harmonic, and for all y G D 

Ny^xnidv) ny{XD e du.Xo ^ 0) = -i„{y)R{dv). 

In addition, Salisbury and Sezer (2012) construct a measurable strictly positive kernel Kn{v] dvi,dv2, . ■ . , dvn) 
from V* to (F*)", concentrated on {(z^i, . . . , Vn) : vi + ■ ■ ■ + = i'}, and an P-nuU set Vq, such that 

^ r ^-{fi,u)+u{xo) ^ if = 

^'^''^"l E,?=i/^e-W>A;,(z.;di.i,...,d^.„)(/i,7.,)(M,7..>---(/^,7.J, if ^ ^ 
is extended X-harmonic for each ^ Vq, is a version of the density in (3) for each p, and also satisfies 

^,{y) = NyiH-'{XD')) 
for every y, every D' D such that y E D' , and every i/ ^ Vq. 
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Proof of Theorem 3. Fix a > 0. It suffices to show the statement of the theorem for /i satisfying (/.j, 1) < a. 
We may assume v ^Q. 

Because D is regular there exists a sequence of smooth domains Di D2 ■ ■ ■ D exhausting D. Let 
xq be the reference point fixed earher, and let fiQ be the point mass at xq. Let fi be another finite measure 
compactly supported in D. Then there exists a smooth domain D G {-Di, D2, ■ ■ ■} containing the point xq 
and the support of fi. Choose H'^{^), Vq, -fi, and Kn as above. 

Let Pi){x, zi, . . . , z„), and Pfj{x, zp), C C {1, . . . , n}, n > 1 be the family of functions defined by equations 
(1) and (2). Let „ £; be the measure on dD defined as in (5). Note that by Theorem 4, 

dp^,,n.b , . _ E^(„)(m,p(-,^Ci))---(/^,p(-,^cJ) 

l^i, . . . , z„j — — = — . {() 

By Theorem 1 there exists A > depending only on D and C U {xq} such that 

p{x,zc,) < A''^''p(a;o,2:cJ for every x £ C. (8) 
Equation (7) and the inequality (8) imply that 



< («A)". (9) 



Let H^{p) — e^^'"^iJ^(yLt). Using the formula (6), the Fubini theorem, and (9), we obtain 

°° f I 

. — 1 J 



n=l 
00 



°° /■ 1 

< E / dvi,..., di^n)Pf,„ „ _D ((aA)"7^i • • • 7i'„) 

/" 1 

= E —^n{v;dvi,...,diyn)Pf_i,„{{aXX^,-f^^)---{aXXf„-f^J) 

.,-1 



We will now show that there exists a Vi C V* such that Vq C Vi, R{Vi) = 0, and for all v ^ Vi, 
P^g{H'^ (kXi).)) < 00 for all integers fc > 1 and i > 1. This will establish the theorem, since for all ly ^ Vi, 

H^{p) = P^(i7"^(X5)) < P^{H^{Xi,)) < P^.iH'^iaXXf,)) < P^„{H'' {kX q)) < 00 

whenever k > aX. Note that H'^{p) > by definition, since the 7^. > 0. 

Fix i and k. Let M > 0. For fixed /i, H'^{p), ^ is a version of [v), so 

Jv 



= P^o Pfcx„, (e- '^^^ Ix.^o)) 

< P^o(e<^^--">Pfcx„^(e-<^°''^^>)) 
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where um is the function on D defined by um{x) = \ogPx{e (^oM)y 

The above derivation shows that Jy. e''^"'^^'^ P^,,XH'' {kXD,))R{dv) < P^o(e<'^°--^("-""'>). This is in fact 
finite if M is sufficiently large. To see this wc use two facts: First, that um converges to u uniformly on 
the closure of D, as they are both continuous functions and um ^ u pointwise in D as M oo. Second, 
it is proven in Salisbury and Verzani (1999) that there exists an e > such that P^aifi'''^'^^''^'') < oo. Hence 
choosing A/ large enough that sup^ u — um < f will ensure that Jy, e~^'^'^^'> Pf^^ {H'^ {kXD^))R{dv) is finite. 

Because e"^"'*^^ is always nonzero, it follows that Png{PP (kXij.)) < oo, for R-a.c. v. Now let Vi be the 
subset of V* defined as 

oo oo 

Vi^Vaij\J\J{veV* : P^.XH'-'ikXD^)) = oo}. 

i=l k=l 

Then = 0, as required. □ 



4. Proof of Theorem 1 

We give the proof first for d > 3, and then describe how to modify it for d = 2. Let 6 ~ dist{C,dD). We 
can find a sequence of smoothly bounded domains Dn with compact closures C„ such that Z) D Co D Ci D 
C2 D . . . D C and 

(n + l)2 

From Harnack's inequality there exists a constant > depending only on Cq and D such that 



dist(C„+i, > ^ for all n>l. (10) 



k{x,z) < (/)fc(xo, 2;)for all x,xq € Cq, and z G dD. (11) 

The 3-G inequality of Cranston, Fabes and Zhao (1988) gives a constant 6 > depending only on D such 
that 

9{x,y) ^ y\'^-'^ + \x - xqI'^-'^ 

— -<0 for all x,Xo,y e £>. (12) 

g{xo,y) 9{x,xo) 

Set (3 = dist{Ca, dD) A 1. Recall that there is a constant Cd such that g{x,y) ~ Cd\x — j/p"** — hx{y) where 
hx is harmonic on D with boundary value Cd\x — Therefore 

g{x, y) <Cd\x- y\'^^'^ for x,y e D, and (13) 
gix, y) >Cd\x^ yp-^ - CdP^-'^ if X, yeD, B{x, /3) C D, (14) 

since Q < < CdP'^^'^ ■ This will imply the existence of finite constants K and _B, depending only on D and 
Co, such that 

1 

< B and (15) 

< _ftr, for all X, xq € Co, x ^ xq. (16) 



gix,xo) 
Cd\x - xoP^'' 
g{x,xo) 



To sec this, note that g{x,xo) is jointly continuous off the diagonal, and everywhere > 0. Thus there is a 
strictly positive lower bound for g on Co x Co less any neighbourhood of the diagonal. (14) gives such a 
bound on a neighbourhood of the diagonal in Cq x Cq. Thus (15) holds. 
To see (16), if x, Xq G Cq and |a; — xqI > § then 

1 ^ 1 „ ,2- 



Scd V2. 

While if X, xq G Cq and |a; — xo| < f then 



g{x, xo)[l + BcdP^-'^] > g{x, xq) + Cdl3^-'^ > cd\x - xof''^ 
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by (14). Therefore (16) holds with K = max(Bcd(/3/2)2-'^, 1 + Bcd^^-''). 

Take N > 1 to be a number whose value will be determined later. Let A = {1, 2, . . . , n}, where n = | A| > 2. 
Let X and xq be two distinct points in C„. Note that 

9{x,y)p{y,ZB)p{y,ZA-B)dy = / g{x,y)p{y, ZB)p{y, ZA-B)dy 

JB(x,-^) 

{x,y)p{y,ZB)p(y,ZA-B)dy. (17) 

We will obtain an estimate for each term in the above decomposition. 

We start with the second term. By the 3-G inequality (12), we have the bound 

' / i\r 2\ '^-2\ 



g{x, y) < g{xo, y) ^^^^ \^\x - xo\' + l^-^ j j bi y e D ~ B{x, ^). 

Combining this with the estimates (15) and (16), wc get 

g[x, y) < g{xo, y)^{— + ^d^^ ) iov y e D - B{x, j^). 

Writing K =^ 6^ and B ^ 6> ^j^!~\ we have 

/ g{x,y)piy,ZB)p{y,ZA-B)dy < [K + Bn'^'-'^^'^n / g{xo,y)p{y,ZB)p{y,ZA-B)dy 

< (a' + ^n2('^-2)^ g^^^^ y^p^y^ ZA-B)dy. (18) 

To treat the first term in the decomposition (17), we introduce a new variable y' = x^ + M{y — x), 
y € B(x, j^)-, for some M > 1 whose value will be determined later. Since y = x + -^{y' — xq), we have 
that dy = M-'^dy', and y e B(x, j^) if and only if y' e B(xo, ^). So by (13) 

/ 9{x,y)p{y,ZB)p{y,ZA-B)dy <cd I \x - y\^~'^p{y,ZB)p{y,ZA-B)dy 

Jb{x,-^) Jb{x,-^) 
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I Iji^dl^o '^piy,ZB)p{y,ZA-B)M '^dy' 



Cd f TH\^o-y'\'^ '^p{y,ZB)p{y,ZA~B)dy'. (19) 

Jb{xo, 



MS , Af2 ' 



Take N = 2M. Note that > 1, so this is consistent with the specification we gave earlier for N. Since 
xo € Cn, we have B{xo, ■^^) = B{xo, ^) C C„_i C Co by (10). Therefore we may apply (16) and (19) to 
get 

/ 9ix,y)p{y,ZB)p{y,ZA-B)dy < g{xa,y')p{y, ZB)p{y, ZA-B)dy' . (20) 

Choose AI > \J2K V 1 and combine the inequalities (20) and (18), to sec that 

g{x,y)piy,ZB)piy,ZA~B)dy < ( K + Bn'^'-'^^^A / g{xQ,y)p{y,ZB)p{y,ZA-B)dy 

< ^ ^ J D 

gi^o, y')p{y, ZB)p{y, ZA-B)dy'. (21) 



B(xo,-^) 



Recall in the above inequality that x and xq are arbitrarily selected points in C„, n > 2, and the constants 
K, B arc independent of x, xq and n. 



To complete the proof we will proceed with an induction argument. Let A be large enough that A > 
(the constant in (11)) and ^> K + Bn'^^'^-^^ for aU n > 1. Then we have that 

p{x, z) < (j>p{xQ, z) < Xp{xo, z) for all x, € Ci, z £ dD. 

Let n > 2, and assume inductively that the following holds for all A; < ?i — 1, x, € Cfc, zi, . . . G dD: 

p[x,zi, ...,Zk)< X''p{xo,zi, . . .,Zk). 

Put A ^ {1,2, ... , n}. Using the inequality (21) for all x, xq in d we have that 

A" f 

g{x,y)p[y,ZB)p{v,ZA-B)dy < — g{xo,y)p{y, ZB)p{y, ZA-B)dy 

D ^ Jd 



9{xQ,y')p{y, ZB)piy, ZA-B)dy'. (22) 



Note that if y' e B{xo, ^) then y = x + ^{y' - xq) is in B{x, ^m^)- Since M > 1, B{x, ^m^) C C„_i by 
(10). Note also that B{xo, C C„_i. By the inductive hypothesis, if y S B{x, jMn^), and y' G B{xo, 
then 

p{y,ZB) < \^^^p{y',ZB) and p{y,ZB) < X^^^^^ p{y' , za-b) 
for all proper subsets B of A. Substituting this in (22) we get that 

A" f 

9{x,y)p[y,ZB)p{y,ZA-B)dz < — g{xo,y)p{y,ZB)p{y,ZA-B)dy 

D JD 



g{xo, 2/')Al^l Al^-^lp(y', ZB)p{y', ZA-B)dy' 



< A" / g{xo,y)p{y,ZB)p{y,ZA-B)dy. 



D 



Hence 



p{x,zi,...,z„) = X! / 9{x,y)p{y,ZB)p{y,ZA-B)dy 

- yi / 9{xo,y)p{y,ZB)p{.y,ZA-B)dy 



B(lA,B^$.A' 

= A'V(xo,zi,...,z„). 

This establishes the result, when d > 3. 

Essentially the same argument works when d = 2. The 3G inequality in that case is classical; We will use 
the version stated in Theorem 1.1 of McConnell (1990), namely that there exists 9 for which 

9ix,y) . ^ gix,y)+g{x,xo) + 1 . ^ ^ „ .„„x 
r- < a r for all x,Xo,y € JJ. (26) 

g[xo,y) 9[x,xo) 

Instead of (13) and (14) we have constants c and C2 such that 

gix,y) < C2log(l/|a:: - y\) + chr x,y e D, and 

g{x, y) > C2 \og{l/\x - y|) - C2 log(l//3) ii x,y E D, B{x, (3) C D. 

As before, this gives (15) and 

C2 log(l/|a:: - XqD 
g{x,xo) 

for aU x,xo G Co, x ^ xq, with K = max(i?C2 log(2/^), 1 + i?C2 log(l//3)). 

9 



< K 



Let > 1 and let x, Xn S C„ be distinct. By (23), we have for y e Z? — B(x, jf^) that 

\ 9{x,xo) J 

< g{xo,y)e(^B{l + c + c2 log(l/|a; - y\)) + l) 

< 5(xo, y)e(B{l + C + C2 \og{Nn^/5)) + l) . 
So for k = 9{b{1 + C + C2 \og{N/5)) + and B = 2eBc2 we get 

/ 9{x,y)p{y,ZB)p{y,ZA-B)dy < {k + Blogn) / g{xo,y)p{y, ZB)piy, ZA~B)dy 

in place of (18). 

Take N = 2M and set y' = xq + M{y — x) as before. In place of (19) we obtain that 



Bix,-^) JB{x 



9{x, y)piy, ZB)p{y, ZA-B)dy < (c2 log(l/|a; - y\) + c) p{y, ZB)p{y, ZA-B)dy 

^{c2 log(A//|xo - y'\) + d^p{y,ZB)p{y,ZA-B)dy' 



^ f K + Bid + C2 \ogM) 

< / 775 g{xo,y )p{y,ZB)p{y,ZA-B)dy . 

Jb{x,,^) 

Now choose M so large that -^'+-^1'=+^'=^^ '°g < 1^ ^j^fj \ go that ^ > k + Blogn for n > 1, and the proof 
proceeds as before. □ 
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